Nuclear modification factor for gluon jets 
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It is shown that the nuclear modification factor can be smaller that unity for jet production at 
small X and at large transverse momentum without any violation of the factorization theorem and 
the initial state effects is able to explain the nuclear modification factor of the order of the one 
measured by ATLAS. In other word, initial state effects are able to describe the jet quenching for 
the gluon jet production. 
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INTRODUCTION 

Recently, the ATLAS collaboration [l| has measured the NMF factor for jet production. It turns out that the 
NMF does not depend on the size of the cone, in which hadrons from the jet decay were measured, and its value is 
considerably smaller than 1 (see Fig. [T|). From Fig. [T] one can see that the suppression is larger for the events with 
small centrality. Such an independence on the size of the come allows us to assume that all hadrons from the jet_decay 
were measured in the ATLAS experiment. If it is so, these da 
discussion below) which is one the most solid result of QCD. 

These data encourage me to ask two theoretical questions. The first question: could we obtain the NMF less that 
unity for production of gluon at very high energy (small x) and with large value of transverse momentum, without 
violation of the factorization theorem? The second one: can be the value of NMF of the order of the value measured 
by ATLAS? Wc would like to emphasize that the goal of this letter to answer these two theoretical questions but not 
to describe the ATLAS experimental data. 
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FIG. 1: Nuclear modification factor(NMF) for jet production. The figures are taken from Ref. [i|] 



Let us first discuss the factorization theorem and why this theorem at first sight leads to the NMF equal to 1. 
Indeed, the factorization theorem states that^the inclusive cross section for gluon production with rapidity Y' and 
transverse momentum px can be written as[! 



<P<j{Y',pt) 

dY'dPpT 



<7ha 



)xiGpr{xi,fiF)x2GTrix2,fJ'F) whcrc Xi 



— =e and X2 



PT _y' 



(1) 



where ^/s = W^ is the energy of colhding particles and ^ denotes all necessary integrations while ^t is the factorization 
scale which should be chosen of the order of pr of produced gluon. Eq. ([1]) has a simple meaning, namely that the 
production of a high pT gluon is proportional to the cross section of the interaction of two partons, with fractions 
of energy xi and X2, multiplied by the probability to find them in the projectile {Gpr) and target (Grr)- All 
interactions in the final state as well as corrections to the initial state (for, example, due to shadowing) lead to small 
small contributions of the order of l/p|.f25'|. 

At short distances (large pt) the gluon structure function for a nucleus is equal to 

xGAix,HT) = AxGp{x,fj,T) (2) 

where Gp is the nucleon gluon structure function. Therefore, 

d^a{A,+A2\Y',pT) ^ d?a{P + P\Y',pT) 

dY'cPpT ' ' dY'cFpT ^^ 

where d^a [P + P\Y' ^pt) dY'd'^pT is the cross section of inclusive production in proton-proton collisions. 
In other words the nuclear modification factor (NMF) Raa, which is defined as 

rfV(Ai+A2|F' ,Pt) 
^AA = — w2„.rP-i-PIV' „^1 = 1 (4) 



N. 



d^a{P+P\Y',pT) 
CO" dY'd-'pT 



where Ncoii is the number of collision that is equal to Ai A2, in the case of large pt- 
The main result of this letter is that Eq. ^ should to be replaced by 

xGA{x,fil/Ql{A;xo)) ^ AxGp {x, fil/Ql{A;xo)) ^ AxGp (^x,fil/(^A^/^Ql{P;xa))^ for ^f « Pt > Qs{A;x) 

. (^) 
where Qs {A;x) {Qs {P;x)) are the saturation momenta ed Y = ln(l/a;) for the nucleus and the proton, respectively. 

Yq = ln(l/a;o) is the rapidity at which the evolution starts. In other word, we state that the linear DGLAP or BFKL 

evolution should be started at different initial condition at low energy (Y = Yq). 

It follows from Eq. (O, that Eq. ^ reads as 

d'^a(^Ai+A2\Y' ,pt) 

_ dV'd'^pT _ a^iGp {xi,fij,/Qi {A; a;io)) X2Gp {x2,fJ.p/Qi {A; xzp)) 

^^ - Ar^^„ if£^;gZ:M - x,Gp{x^,^ll/QUP;x^o))x^Gp{x2,^l%/QUP■,^2o)) ^' 

For the sake of numerical estimates we will derive the QCD factorization formula from the fcy-factorization, in 
order to specify the factorization scale fip and all numerical coefficients in Eq. ^ . Using kp factorization [7H10| , the 
inclusive cross section has the form 

d^ ^1^1 '^'^^ '^^' (^^' ^^0 ^^' (^^'^^^ " ^0 ^^^ 

where (/)q* is the probability to find a gluon in the nucleus {Ai), that carries the fraction of energy xt with kp transverse 
momentum, and ols = UsNc/tt, where Nc is equal to the number of colours in the SUn^ colour group. At large values 
of pp we can re- write Eq. ([7]) in the form 

da 2™, r r^ ,,, ,4, / r \ ,4„ , . . ,4, , . . r^ 



dy d^pp p^ 



d^kpcj)^^ ixi]kp\ (f>^^ {x2]Pp) + 4>G^ {xi]pp) / d^kpcf)-^" (x2;kp] | 



27ra, d 



2 ■^[^iGai{xi,Pt)x2Ga2{x2,Pp)J (8) 

For the purpose of deriving Eq. (|S]) , we used the relation between the un- integrated gluon density (f> and the gluon 
structure function (see for example Ref. [13i]): 

xG{x,Q^) = / (t){kl)dkl (9) 



Using Eq. (|5]) one can see that Eq. ^ can be rewritten as 



^(xiGp {xi,pl/Ql (A;xio)) X2Gp {x2,p^/Ql {A;x2o)) ) 
Raa = -^ f (10) 



^ Ix.Gp {x^,pI/QI (P; a;io)) x^Gp {x^^vl/Ql (P; x^^)) 



RELATION a;GA(a3,p|/Q? (A; xo)) = AxGp [x,p'^/ [A'-'^Ql {P;xo) 



The proof of Eq. ([Sj is based on two observations. First, the equations in the region of low x can be written as 
an evolution in rapidity, both in the approach of DGLAP ll| and in the BFKL 12] approach. They both have the 
following general form 



dxG{Y;pD 
dY 



/ d?kT K {pT, kr) xG (F, k^) (11) 



where Y = ln(l/a;). 

Second, since for large pt we can restrict ourselves to the linear evolution: DGLAP or BFKL, one can see that the 
solution of Eq. ([TT|) can be written as 

xG{Y,^=\n{pl/KlcD)) = j di'G{Y-Y^,i-C)xG,^{Y^Y^,0 (12) 

where G {Y — Y^,^ — ^') is the Green function of Eq. (|lip : the solution to Eq. (|lip with the initial condition 
G (Y ^Y{),S, — ^') = S {^ — £,'). xGin (Y — Yoj C') is the initial condition for the gluon structure function. 
For linear evolution the solution to the equation reveals the following property of the propagator 13l4l5l| : 

xG(y, i = In {pI/KIcd) ) = I d^' G{Y^Y',^^ ^ xG {Y' , O (13) 

for any value of Y'. We can prove Eq. (fT3)) by solving Eq. (TTTI) using the double Mellin transform 

^G(r,0 = / ^ ^e"^+-«5(^,7) (14) 



£ — joo ^""^ Je~ioo ^^* 



Plugging Eq. (fT4|) into Eq. ([TT|) we obtain the solution 

LU — UJ (7) (15) 

where w (7) is the Mellin image of the kernel in Eq. (fTTj) . It can be verified that the Green function has the form 

GiY-Y',0 - / ^e'^(^)(^-^')+7(C-«') (16) 

Je-ioo 27ri 

while the solution of Eq. (J12p can be written in two equivalent forms 

xG{Y,0 = P";^e-(-)^+'^«5»(7) = P°°|^e-^+-^(-^)«5™(c.) (17) 



where g^ is the Mellin image of xGin {Y — YqjS^') in Eq. ()12p and 7 (w) is the anomalous dimension which can be 
obtained as solution of Eq. (fT5|) . 

Using Eq. p^ and Eq. p7)) for the r.h.s. of Eq. ([T^ . and integrating over ^' we obtain the l.h.s. of this equation. 

Our standard approach to low x evolution consists of two steps. First, we assume that there exists a small x = xq 
(Y — Yq — ln(l/xo)) which is large (Fo ^ 1), but at the same time small enough such that cXsYq ^ 1 . For this type 
of Yq we have the theoretical formula (the McLerram-Venugopalan formula [16] ) for the scattering amplitude of the 
dipole, which looks as follows in the most simplified form 

Nir,b;Yo) = 1 - c^p {-r^Ql{Yo,b) /4) (18) 



where N is the imaginary part of the scattering amphtude; Q^ is the saturation scale atY = Yq , r is the dipole size and 
b is the im pac t parameter. The main features of Eq. (|18p are that the amphtude of Eq. p8|) has the geometric scahng 



behaviour [l7| at low energy {Y = Yq), since it is a function of only one argument N (r, b; Yq) = N {t = r'^Qs {Yq, b)) . 

The second step is to use the evolution in the region Y > Yq. In our case, when the typical size of the dipole is 
small (r w 1/pt), we can use the linear BFKL or DGLAP equations. It should be stressed that we use the rigorous 
theoretical formula of Eq. (|18l) , which shows the geometric scale behaviour only for the initial condition, for the linear 
evolution. 

We wiU show below that for the initial condition xqG {Yq, kx) = S Ql{YQ)g (r = k'^j./Ql {Y)). Using Eq. dTS]) with 
Y' = Yq we have 

xG(y, e = In (p't/AIcd) ) - J d^' G{Y-Y,,C- O xqG (Fq, C) 

= SQl {Yo) J ^G {Y - ro,p|/4) g {k^/Ql {Y,)) 

where r = k'^/Ql and S is the area of the target. In Eq. ([T^ we used Eq. (IT51) . The function 7(7) does not depend 
on the character of the target, and it is the same for the nucleus and for the nucleon if we consider Eq. (|18l) for both. 
Since Qs (A; Fo) = A^l^ Qs (P; Yq) and Sa = A^/^ , one can see that Eq. © follows from Eq. (fTg]) . 

Now let us show that Xf^G {Yq^Ut) — S Ql(YQ)g{T). The relation between (pQ ix^kT} and the dipole scattering 
amplitude N {x, r; b) has the form [l^ 

0^(xo;fcT) - -^^ j (fb(fTe'^^-^^lN^{Yo;r-b) (20) 



where 



N^{Y^-r,b) = 2N{Yo;r,b) - n' (Yo;r,b) ^^^ 1 - exp (-r^g^ (Fq; 6) /2) (21) 



Plugging Eq. (1^ into Eq. (^01 and using the simplest model for the b dependence, viz. Ql{Yo;b) — 
Qs (Yo) Q {Ra — b) we obtain 

c^^ixo-M) ^c^^{T = kl/Ql{Ya)) = 2^5rexp(-2T) (22) 

Using Eq. ^ we calculate the gluon structure function at y == lo- It is equal to 

xoG(xo,p|) = \SQ^{Y)) (1 - (l + 2T)cxp(-2r)) = SQl{Yo)g{T) (23) 

Therefore, we arrive at the result xqG {Yq, kr) = S Ql{YQ)g (r) with the particular function g (r). The corresponding 
function / (r) (see Eq. ([T^) turns out to be equal to / (r) = 2T (2 — 7) /j. 

Eq. p9| is the main result of the paper. It should be stressed that solution which we found, does not lead to the 
geometric scaling behaviour since it depends on both variable: x and pT- However, Eq. (fT9|) shows that pT in the 
DGLAP evolution scales with different momentum that depends on the target. In this sense this equation shows 
the PT scaling behaviour. The fact that in DGLAP evolution px is scaled by the typical momentum in the initial 
condition is well known and follows directly from the conformal symmetry of the DGLAP equation at fixed QCD 
coupling. Two ideas are new in the proof of Eq. ^T9\\ : the evolution at low x has to start with rather low initial xq 
and the initial condition at this x — Xq has to show the geometric scaling behaviour. Both of these results follows 
from the Color Glass Condensate/saturation approach and are rigorous results for QCD. 

To illustrate the situation we consider the semi-classical approach to the solution of the DGLAP equations [l9J- In 
this approach we take the integral over 7 in Eq. (|19p using the steepest decent method and the solution is characterized 



Y=ln(1/x) 

Q=Qs{p,x) Q=Qs(A,x) 
CQG/saturation 



HV 




ln(Q?/^cD) 



FIG. 2: The semi-classical trajectories for the nucleus (solid lines) and nucleon (dashed lines) gluon structure function in the 
region of large and small x. In the Color Glass Condensate (CGC)/saturation domain the structure function has a geometric 
scaling behaviour. 

by the trajectory: the hne in (ln(l/a;),ln (Q^/^qcd) pl^-ne. The trajectory shows what values of xq and Qo in the 
initial condition are essential to find the solution in the point (ln(l/.T),ln {Q^ /^qcd))- Fig- [H shows the trajectories 
for gluon structure function in the region of low x and in the region of x ~ 1. These trajectories are denoted by solid 
and dashed lines for nucleus and nucleon targets, respectively. On can see that for a; ~ 1 both trajectory for nucleus 
and for proton targets started at the value oi Q = Qo which corresponds to the same scaling of pr for interaction with 
nucleus and nucleon. However, for the region of low x the situation changes crucially and the evolution starts from 
the value of Qo that should be determined by the initial condition given by the Color GlassCondensate/saturation 
approach. Fig. [2] illustrates the main result of the CGC /saturation approach: we cannot use the initial condition at 
fixed arbitrary Q = Qq since its contradicts the unitarity constraints. 

ESTIMATES FOR NMF 



At first sight the value for Raa could be calculated in a direct way, just by using Eq. ^ and various different fits to 
the available DIS data, using the DGLAP evolution equations. By plugging in Eq. ^ the gluon structure functions 
that are available on the market (see Ref. 21 1) and using Q^ {A; Yq) — A^^^Ql {P; Yq) , then we obtain Raa for the 
gluon jet production at the LHC, in the central rapidity region for the lead-lead collisions. This estimates are shown 
by the black line in Fig. [31 In Fig. [3] are plotted the NMF using the MRSTW parameterization [20| but we check that 
the other parameterizations give approximately the same value for Raa- 
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FIG. 3: Nuclear modification factor (NMF) for jet production. The black line shows Raa from Eq. ^ for MRSTW param- 
eterizations of the gluon structure functions (see Ref.[2y]. Raa from Eq. (|25|l is plotted by the blue line, cis is taken to be 
as = 0.2 and a;i = a;2 = 10"^ 



However, this type of approach suffers from two major defects. First, in all parameterizations the running QCG 
coupling is used. The transverse momenta pr enters to the expression for the running QCD coupling as In {pr/ ^qcd) 
and, therefore, they are not rescaled with the saturation momenta. Second, the gluon structure functions, that we 
used, where extracted from the experimental data using the DGLAP evolution in transverse momenta, but not in Y . 
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FIG. 4: Nuclear modification factor (NMF) for jet production. The data is taken from Ref.[lJ. The curves are calculated from 
Eq. (I26[l while the value of the saturation momenta for each centrality events are estimated in the KLN- procedure 22], with 
the number of participants taken from Ref.[l|. The peripheral collisions are the events with centrality 80 -i- 100%. The values 
of xi and X2 are chosen to be xi = X2 = 10~^. 



We believe that we can obtain more reliable estimates by using Eq. ([T^ . Eq. P7)) and Eq. (1^51) . It is known that the 
leading order DGLAP anomalous dimension can be approximated with an accuracy of 5% by the following expression 

M 



7 = "sf Ij or a; (7) 



7 + "s 



(24) 



Substituting a; (7) from Eq. (|24p and taking the integral using the steepest decent method, we obtain the following 
answer: 



xG {Y, T ^) = f (7sp) exp f 2 ^/a^ Y In (r) - a^ In (r) j with jsp 



I a.Y 
ln(r) 



(25) 



where r = Px/Ql (^0) and Y = ln(l/a;). 

The blue curve in Fig. [3] describes Raa from Eq. (|10p with the gluon structure function given by Eq. (|25|) . This 
estimate leads to smaller values of Raa = 0.4 ~ 0.5 (see Fig. [3]). 

Using the KLN approach [22|, we can estimate Rcp which is defined as 



Rcp 



dN" 



^j^periph 



jycentr ^YiPpT I N^Z"^^ dYcPpT 
'^^i''' coll 



(26) 



dp'^ 



[xiGp {xi,p'^/Ql {centr;x 10)) X2Gp {x2,PtIQI {centr; x 20)) 



■j^ VxiGp {xi.p^/Ql {periph; xw)) X2Gp {x2,p\IQ'\ {periph; X20)) 

where N^^Pi*^{Q {centr ;Yo)) and N^^ii^ {Q {periph; Yq)) are the number of collisions ( saturation momentum) in 
the events with fixed centrality and for peripheral collision, respectively. Based on the KLN approach we know the 
saturation momentum for the events with fixed centrality, and using Eq. (PB|) we calculate the Rcp (see Fig. S]). 



CONCLUSIONS 



The main result of the paper is that the two questions which have been formulated in the introduction, namely, 
could we obtain the NMF less that unity for production of gluon at very high energy (small x) and with large value 
of transverse momentum, without violation of the factorization theorem; and can be the value of NMF of the order 
of the value measured by ATLAS, have the affirmative answers. 

Eq. ([5]) shows that the NMF for the gluon jet production turns out to be less than 1. This equation is proven using 
two main ingredients from the Colour Glass Condensate (CGC)/saturation approach; the typical xq in the initial 



condition for the evolution in the region of low x even at high values of pt is rather low and we need to use the 
Color Glass Condensate approach to describe them; and the initial condition at fixed Y = Yq depends on one variable 
r = Pt/Q^O^o), and this initial condition has the same form for the scattering with nuclei and protons. The last 
assumption is essential for the estimates of Raa, but not for Rcp- For Rcp we can use a weaker assumption, namely, 
that the McLerran-Venugopalan formula can be used for the description of the peripheral collisions. 

As have been claimed in the introduction we do not pretend that we can describe the data since we consider the 
theoretical example of very small and fixed x and large values of pt ■ The kinematical region were the experimental 
data are taken, perhaps, is quite different. Our estimates show that the scale of the effect is of the oder of the 
experimental one. Therefore, the fact that we do not reach a good agreement with the experimental data does not 
discourage us. It should be mentioned that the simple formula that we use (see Eq. (|25p ) relies on the saddle point 
approximation, and on the particularly simple form of Eq. (|18|) for the McLerran-Venugopalan formula. 



For serious comparison with the experimental data we have to develop an approach similar to one in Ref. j24j . We 
believe that we find the explanation, why NMF calculated in this paper as well is in many others based on CGC 
approach , lead to the value of the NMF less than unity. 

The conclusion from the paper can be formulated in one sentence: the NMF can be smaller that unity for jet 
production at low and fixed x and at large transverse momentum pr without any violation of the factorization 
theorem and the initial state effects are able to explain the NMF of the order of the one measured by ATLAS. In 
other word, initial state effects are able to describe the jet quenching for the gluon jet production. 

This mechanism leads to stronger suppression at smaller value of pt (see Fig. [3]-b) and it should be taken into 
account in the explanation of the NMF for produced hadrons. 
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